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Abstract. The beaming angle ( m is the main uncertainty 
in gamma ray burst energy requirements today. We sum- 
marize predictions for the light curves of beamed bursts, 

and model the R band light curve of GRB 970508 to de- 



rive ( m 30 . This yields an irreducible minimum energy 
requirement of 3.4 x 10 49 ergs to power the afterglow alone. 



1. Beamed Gamma Ray Burst Afterglow Models 

In the Rome meeting I presented a derivation of the dy- 
namical behavior of a beamed gamma ray burst (GRB) 
remnant and its consequences for the afterglow light curve. 
(Cf. Rhoads 1999 [Paper I]). Here, I summarize these re- 
sults and apply them to test the range of beaming angles 
permitted by the optical light curve of GRB 970508. 

Suppose that ejecta from a GRB are emitted with ini- 
tial Lorentz factor Tq into a cone of opening half-angle 
Cm and expand into an ambient medium of uniform mass 
density p with negligible radiative energy losses. Let the 
initial kinetic energy and rest mass of the ejecta be Eq 
and Mo, and the swept-up mass and internal energy of 
the expanding blast wave be Af acc and Ei nt . Then energy 
conservation implies TE[ nt w r 2 M acc c 2 w Eq w constant 
so long as l/r £ Af acc /M <, T . 

The swept-up mass is determined by the working sur- 
face area: dM acc /dr » 7r(C m r + c s t co ) 2 , where c s and 
t co are the sound speed and time since the burst in 
the frame of the blast wave + accreted material. Once 
r 1/Cmj c s t co > £ m r and the dynamical evolution with 
radius r changes from F cx r~ 3 / 2 to T a exp(— r/r r ) 
(Rhoads 1998, & Paper I). The relation between ob- 
server frame time and radius r also changes, from 
t® ex r 1 / 4 to t® oc exp(r/[2r r ]). Thus, at early times 



r oc i 



-3/8 



while at late times V cx i 
,2 „\V3 



1/2 



The charac- 



teristic length scale is r r = (Eo/ircjpY'^ , and the charac- 
teristic observed transition time between the two regimes 

is t e , b « 1.125(1 + z) (E c 3 /[pc 8 s C}) 1/3 Cm 3 , where z is 
the burst's redshift. 

We assume that swept-up electrons are injected with 
a power law energy distribution N(£) cx £~ p for £ = 



7 e m e c 2 > Exam ~ Ce^pC 2 r, with p > 2, and contain a 
fraction £ e of £"i n t ■ This power law extends up to the cool- 
ing break, £ coo i, at which energy the cooling time is com- 
parable to the dynamical expansion time of the remnant. 
Above £ coo i, the balance between electron injection (with 
iVj n j cx £~ p ) and cooling gives N(£ ) cx £~( p+1 \ 

We also assume a tangled magnetic field containing a 
fraction £ B of E- lnt . The comoving volume V co and burster- 
frame volume V are related by V co rs V/T cx M acc /r, so 
that B 2 = 8Tr£ B E int /V co cx T 2 and BcxT. 

The resulting spectrum has peak flux density F v ,®, m oc 
TBM acc / max(C m ,r -2 ) at an observed frequency ^©, m cx 
TBS mln /(l + z) cx T 4 /(l + z). Additional spectral features 
occur at the frequencies of optically thick synchrotron 
self absorption (which we shall neglect) and the cooling 



frequency 



(which is important for optical obser- 



vations of GRB 970508). The cooling break frequency 
follows from the relations 7 coo i ~ (67rm e c)/ '(cr T TB' 2 t S) ) 
(Sari, Piran, & Narayan 1998; Wijers & Galama 1998) and 
^©,cooi oc Ti?7 2 ool oc (r 4 *!) -1 . I n tne power law regime, 



Fv,<s, m oc i®, v®, m oc t B 3 /2 , and z/, 

the exponential regime, F v< ® im cx t^ L , ^®. m oc and 
^©,cooi oc tgg. The spectrum is approximated by a bro- 
ken power law, F v cx v~P , with fi w —1/3 for v < z/® im , 
P w (p - l)/2 for ^ ffi!m < j/ < i/© )0OO i, and /? w p/2 for 

V > ^©,cool- 

The afterglow light curve follows from the spectral 
shape and the time behavior of the break frequencies. 
Asymptotic slopes are given in table 1. For the T ~ 1/Cm 
regime, we study the evolution of break frequencies nu- 



i,cool OC t s 
-1 



-1/2. 



while in 

2 



merically. The results for v§ m and F v 



are given m 



Paper I. For i/g 



.cool J 



a good approximation is 



V®, cool — 

1 

X 



5.89 x 10 13 (W*e,fc) 1/2 + i- 34 x 1Ql 



Hz 



I : / V c/a/37 

3 \ -5/6 



17/6 / ^ N -3/2 

SB 

0.1 



p ■ cm 



^ /10 53 erg\ 2/ VCn 



a/4 



0.1 



-4/3 



2 
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i/abs < v < u m v m < v < V coo \ v coo \ < v 
te < *e,6 1/2 3/4 - 3p/4 1/2 - 3p/4 

*e > *e,6 -I/ 3 -P ~P 

Table 1. Light curve exponents a as a function of fre- 
quency and time. Here F w ,© oc iX. 



Model 
1 
2 
3 



25.74 
25.55 
25.36 



P 

2.36 
2.20 
2.04 



lg(«e 

9 

5 

3.75 



i/ da y) 



20.32 
20.33 
20.32 



X 2 /d.o.f. 
4.34 
3.80 
3.55 



Table 2. Fitted break times i© ; 6 and magnitudes R c (to) 
(at fiducial observed time to = 3.231 day) for beamed 
GRB afterglow models for three pairs of acceptable host 
galaxy magnitude Rh and electron power law index p. The 
fit included all 43 data points with 1.3 day < t© < 95 day 
in the compilation by Garcia et al (1998). 



2. Application to GRB 970508 

In the best-sampled GRB afterglow light curve yet avail- 
able (the GRB 970508 R band data), the optical spectrum 
changed slope at t© ~ 1.4 day, suggesting the passage of 
the cooling break through the optical band (Galama et 
al 1998). We explore the range of acceptable beaming an- 
gles for this burst by fitting the afterglow light curve for 
1.3 day < i© < 95 day assuming that ^©, C ooi < c/0.7fim. 

The range of acceptable energy distribution slopes p for 
swept-up electrons is taken from the optical colors. Precise 
measurements for 2 day ^ t© & 5 day give F v oc with 
(3 = 1.10 ±0.08 (Zharikov, Sokolov, & Baryshev 1998), so 
that p = 2.20±0.16. We take this value to hold throughout 
the range 1.3 day < i© < 95 day, thus assuming that p 
does not change as the afterglow evolves. We subtract the 
host galaxy flux {R H = 25.55 ± 0.19; Zharikov et al 1998) 
from all data points before fitting. 

We fixed values of Rh and p, and then executed a grid 
search on the break time t©^ and normalization of the 
model light curve. Results are summarized in table 2 and 
figure 1. The final \ 2 P er degree of freedom is ~ 4. 

These large % 2 values make meaningful error estimates 
on parameters difficult. Let us suppose \ 2 is large because 
details omitted from the models (clumps in the ambient 
medium or blast wave instabilities) affect the light curve, 
and so attach an uncertainty of 0.1 mag to each predicted 
flux. Adding this in quadrature to observational uncer- 
tainties when computing % 2 , we obtain % 2 /d.o.f. ~ 1. 
Error estimates based on changes in \ 2 then rule out 
lg(te,b/ day) < 3.5 at about the 90% confidence level even 
for our "maximum beaming" case (p = 2.04, Rh — 25.36). 

To convert a supposed break time i©.b into a beam- 
ing angle Cm, we need estimates of the burst energy per 
steradian and the ambient density. Wijers & Galama 
(1998) infer E /VL = 3.7 x 10 52 erg/(47r Sr) and p = 
5.8 x 10 _26 g/cm 3 . Combining these values with i©^ <; 



log(t / days) 
1 1.5 




Fig. 1. Upper panel: The Cousins R band light curve for 
GRB 970508 with the three fits shown in table 1. Lower 
panel: Residuals for the data and for models 2 and 3 (in 
order of increasing curvature) relative to model 1. A host 
galaxy flux corresponding to Rh — 25.55 has been sub- 
tracted from all data points. 

10 3 5 day gives Cm ^ 0.5 rad w 30deg. Eq/VL and p are 
substantially uncertain, but because Cm °c (p/Eq) 1 / 8 , the 
error budget for Cm is dominated by uncertainties in p 
rather than in Eq or p. 

This beaming limit implies £1 > 0.75 Sr, which is 6% 
of the sky. GRB 970508 was at z > 0.835 (Metzger et 
al 1997). We then find gamma ray energy E^ — 2.8 x 
10 50 erg x (O/0.75 Sr)(d L /4.82 Gpc) 2 (1.835/[1 + z}). If the 
afterglow is primarily powered by different ejecta from the 
initial GRB, as when a "slow" wind (ro ~ 10) dominates 
the ejecta energy, then our beaming limit applies only to 
the afterglow emission. The optical fiuence implies i? op t — 
3.4 x 10 49 erg x (ft/0.75 Sr)(d L /4.82 Gpc) 2 (1.835/[1 + z}). 
The irreducible minimum energy is thus 3.4 x 10 49 erg, 
using the smallest possible redshift and beaming angle. 
We have reduced the beaming uncertainty, from the factor 
~ T 2 ~ 300 2 ~ 10 5 allowed by 7-ray observations alone to 
a factor (4ir Sr)/(0.75 Sr) ~ 20, and thus obtain the most 
rigorous lower limit on GRB energy requirements yet. 
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